Heisenberg magnet with modulated exchange. 
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A modification of the ground state of the classical-spin Heisenberg Hamiltonian in the presence 
of a weak superstructural distortion of an otherwise Bravais lattice is examined. It is shown that a 
slight modulation of the crystal lattice with wavevector Q c results in a corresponding modulation of 
the exchange interaction which, in the leading order, is parametrized by no more than two constants 
per bond, and perturbs the spin Hamiltonian by adding the "Umklapp" terms ~ ^Q^q±Q ■ As 
a result, for a general spin-spiral ground state of the non-perturbed exchange Hamiltonian, an 
incommensurate shift of the propagation vector, Q, and additional new magnetic Bragg peaks, at 
Q ± nQ c , n = 1, 2, appear, and its energy is lowered as it adapts to the exchange modulation. 
Consequently, the lattice distortion may open a region of stability of the incommensurate spiral 
phase which otherwise does not win the competition with the collinear Neel state. Such is the 
case for the frustrated square-lattice antiferromagnet. In addition, the "Umklapp" terms provide 
a commensuration mechanism, which may lock the spin structure to the lattice modulation vector 
Q c , if there is sufficient easy-axis anisotropy, or a magnetic field in an easy plane. 

PACS numbers: 75.10.-b 75.25.+Z 75.50.-y 



I. INTRODUCTION 



An interplay between the distortion of a crystal lat- 
tice and the magnetic properties of the material has 
recently become a subject of renewed interest. One 
problem which provides strong motivation for study- 
ing the effect of a weak superstructural modulation 
on the spin system is that of striped phasesii^ These 
charge-ordered states are found in lightly doped high- 
T c cuprates La 2 _ x Sr :r Cu04 +y (LSCO) and in related 
nickelates, and are always associated with a weak su- 
perstructural distortion of the original "stacked square 
lattice" structure of the un-doped parent material. In- 
commensurate magnetism in these compounds is usu- 
ally interpreted in terms of the segregation of the doped 
charges into lines which separate the antiferromagnetic 
domains ("stripes"), characteristic of the un-doped ma- 
terial. There is also a modulation of the crystal struc- 
ture induced by the charge-stripe segregation, but it is 
often too small to be observed in experiment^ It is clear 
that the essential effect of the stripe order on the spin 
system of cuprates is that of a periodic modulation of 
the exchange coupling in the Heisenberg spin Hamilto- 
nian which describes their magnetic properties: 3 How- 
ever, only the simplest "average" consequence of stripe 
superstructure, in the form of the effective weakening of 
the exchange coupling in the direction perpendicular to 
stripes, has been considered so far^ A similar problem, 
of an interplay between the spin order and a coopera- 
tive Jahn- Teller distortion accompanying the charge or- 
der, arises in the context of the charge-ordered phases in 
doped manganites^ 

A number of examples not related to charge ordering 
but no less interesting, involve an intriguing interplay be- 
tween small superlattice modulation and spin structure 
in spin-frustrated antiferromagnetic dielectrics. In the 
square-lattice antiferromagnet, a distortion could actu- 



ally be the source of the frustration. For example, it may 
generate a well-known generalized Villain model. 6 In the 
triangular-lattice antiferromagnet (TLA) , a weak distor- 
tion may partially release frustration, and result in com- 
plicated spiral phases. Among the simple realizations of 
the distorted TLA are the so-called "row models" which 
were extensively studied in the past£&2ii£ However, pre- 
vious studies were mainly restricted to a specialized anal- 
ysis of a few particular models; no general approach that 
would allow a unified treatment of the effect of a small 
lattice distortion on a spin system has been developed so 
far. 

A number of interesting experimental examples of 
distorted triangular lattice antiferromagnets (DTLA), 
which instigated this study, are found among the 
CsNiCl3-type compounds with the general chemical for- 
mula ABX3. In the (anti)ferroelectric phases that are 
realized in some of these materials at low temperatures 
relevant for the magnetic order, the P&s/mmc hexago- 
nal symmetry in which they crystallize is lowered, and a 
fully frustrated triangular lattice inherent in the original 
CsNiCl3-type "stacked triangular lattice" crystal struc- 
ture is slightly distorted. Typically, the high-temperature 
hexagonal structure with a Bravais lattice of equiva- 
lent magnetic B 2+ sites is changed to either a hexag- 
onal -P63 / cm structure with a 3 times larger unit cell, 
or to a large-cell orthorhombic structure. These super- 
structures are characterized by the appearance of the 
superlattice Bragg reflections at Q c = (77, 77, 0) fii with 
r) = 1/3, or i] — 1/4,1/8,..., respectively. In some 
cases, as in KN1CI3, both phases are found to coex- 
ist at low temperature^ Perhaps, the most intriguing 
is the case of RbMnBr3, in which most experiments 
find the orthorhombic low-T phase ji&iiiia and an in- 
commensurate spiral spin structure with propagation 
vector Q = (1/3 + q, 1/3 + q, l) , 13 i 14 i 15 i 16 in place of 
the commensurate "triangular" antiferromagnetic order 
with Q = (1/3,1/3,1), which is characteristic of the 
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non-distorted hexagonal materials CsMnBr3, CsNiChj, 
ete i 17 i 18 i 19 In a magnetic field of about 3 T applied in 
the easy plane the spin structure becomes commensu- 
rate, with Q = (1/8, 1/8, 1). In the similar orthorhombic 
modification of KMCI3, which is a related material but 
with an easy-axis spin anisotropy, this latter structure is 
realized already at H = 0. 

Until now, these experimental findings remained to a 
large extent unexplained. One reason for this is that tra- 
ditionally, the effect of each particular lattice distortion 
on the spin Hamiltonian was considered separately, by 
devising a specific, generally multi-sublattice spin model 
(eg, the "row models"), where the distortion simply de- 
fines the particular setup of the near-neighbor exchange 
interactions. For long-period structural modulations, 
this approach leads to models with a large number of 
inequivalent spin sites (eg up to 8 for Q c = (1/8, 1/8, 0)), 
resulting in tremendously complicated spin Hamiltoni- 
ans, and, therefore, the analysis has never been carried 
through. The same problem is outstanding for the stripe 
phases in LSCO cuprates, where the most stable super- 
structure has a pitch of about 1/8 (curiously, it is the 
same as that of the antiferroelectric lattice distortion in 
RbMnBrg— ). In addition, the modulation has an even 
longer period at small doping, and, in general, can also 
be incommensurate. Here we devise an alternative ap- 
proach, which lays grounds for the consistent and general 
explanation of spin incommensurability, commcnsuration 
transition, and other phenomena arising from the lattice 
distortion, that were mentioned above. We treat the ef- 
fect of an arbitrary but small lattice distortion on the 
microscopic spin Hamiltonian in the perturbation frame- 
work. The analysis in this paper most directly applies in 
the case of the dielectrics with localized spins, although 
we expect it to hold also for the doped perovskites, to 
the extent that the itineracy effects can be neglected. 

Consider a system of N equivalent spins on a simple 
Bravais lattice, coupled by Heisenberg exchange interac- 
tions. The model Hamiltonian, which allows also for a 
uniaxial spin anisotropy and a Zeeman energy, is 

h = J2 Ja ^ -s 3 + dJ2 (si) 2 -I> S h (1) 

i,j i i 

where Jy = Jji parameterize the exchange coupling be- 
tween the spins at lattice sites i and j, D is the anisotropy 
constant, and h — g^sH is the magnetic field. With- 
out the anisotropy and magnetic field, the classical 
ground state of P is a planar spin spiralj2fli2ii22i22i2i 
Sj = (Scosdj,Ssindj,0), 9j = (Q ■ rj). The order- 
ing wave vector Q corresponds to the minimum of the 
lattice Fourier transform of the exchange interaction, 
Jq = J2r i3 exp(-ig • Tij), r i3 = Tj - r t (this in- 
cludes ferro- and antiferromagnetism, corresponding to 
Q = and, eg, Q — (tt, tt, it), respectively). For non- 
zero D and h, the spin structure is, in general, modi- 
fied by the appearance of higher-order harmonics in the 
lattice Fourier transform of the spin distribution Sq = 



J2r Sjevp(~iq ■ Tj), with the wavevectors Q n = nQ, 
in addition to Sq?^^^ For small D and h the corre- 
sponding terms in the Hamiltonian Q can be treated 
as a perturbation. The perturbation expansion for the 
correction to the spiral winding angle has the form 69 j = 
( a « cos ( n Q ' r j) + Pn sm ( n Q ' r j)) ' where the coef- 
ficients a n and (3 n are of the order 0((^)^ , (7)™) — 

In this paper we consider what happens to the classical 
spiral ground state of the Hamiltonian if a weak su- 
perstructural modulation of the original Bravais lattice 
appears in the crystal whose spin system it describes. 
It is clear that results obtained for classical spins are 
subject to corrections from quantum and thermal fluc- 
tuations, and these are often crucial. Nevertheless, if we 
find the spiral state that has the lowest energy already on 
the classical level, inclusion of the 1/S corrections may 
still result in this state being the ground state (clearly, 
this will always be the case for large enough spin S). In 
fact, this may appear to be the case even for small spins, 
so long as the spin- wave theory holds. While account for 
the fluctuations is important, it falls beyond the scope of 
this paper and is deferred to further studies. 

II. SUPERLATTICE DISTORTION AND THE 
SPIN HAMILTONIAN 

Consider a slight distortion of the crystal structure 
which is characterized by the appearance of additional, 
weak supperlattice Bragg reflections at wavevectors ±Q C 
in the Brillouin zone of the non-distorted Bravais lattice. 
Most generally, such a superstructure corresponds to a 
small harmonic modulation of the positions, r^ 1 , of the 
ions and the characteristic symmetry points of the local 
electron density distribution (orbitals) in the lattice, 

(r»Y = r» + ^,e» = e» cos(Q c • r£) + e£ sin(Q c • r£). 

(2) 

Here j numbers the sites of the original Bravais lattice of 
magnetic ionsj2& and fi indexes positions of ligands and 
symmetry points of the magnetic orbitals within the unit 
cell. It is shown in Appendix lAl that a leading correction 
to the exchange coupling in this case is most generally 
expressed as, 

J'ij = Jij + hj cos (<3c ' Rij) - fij sin(Q c • Rij), (3) 

where Rij — \(Ri + Rj) is the middle of the ij bond, 
j'^ = j,., j'-j — are real and symmetric, and, most 
importantly, satisfy all symmetries of the original lattice 
which leave Q c and the polarizations e^, invariant; the 
same is true for J^. In particular, Jy , j\- and are 
invariant under the translation group of the non-distorted 
lattice. Finally, Q c = nQ c , where n = 1,2, is the order 
of the leading correction to Jy of Eq. © in terms of the 

small parameter e ~ < 1 U'i^hj ~ <"-A, ; 2 '' 

In the first order in e, lattice modulation results 
in the modulation of the exchange coupling J[, with the 
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same wavevector Q c . It is described by Eq. JSJl, with 
Q c = Q c , Ji 3 = Jij, and jij,j'ij, defined by Eqs. HA16f> . 
(|A17|) of Appendix ^ The structure of the expressions 
for jy is quite illuminating. There are two contribu- 
tions, one of which is ~ sin { ^"2 3 ) > an( l, therefore, an- 
tisymmetric with respect to Q c — > — Q c , and the other de- 
pends on the relative alignment of Q c and the distortion 
polarizations with respect to the bond geometry. Only 
the first contribution survives in the simplest case when 
exchange depends on the bond length alone, = J(ry); 
in this case the bonds that are perpendicular to Q c are 
not changed by the lattice distortion. 

In many important cases the first-order corrections 
vanish, and the leading correction in Eq. @ is ~ e 2 Jij, 
in which case Q c — 2Q C , and Jij,jij are given by Eqs. 
(|A9I) - (|A12|) . In general, the latter have to be amended 
in accordance with l|A14(l . (|A15(I . Importantly, a Q c - 
independent contribution which determines the leading 
correction to the bond strength Jy , 

Jij Jij ~\~ & Jij , 3 Jij ^ ^ , (4) 

always appears in this order. Because the correction to 
the bond strength, 5 Jij, arises from replacing 2cos 2 (Q c • 
Rij), 2 sin 2 (Q c ■ Rij ) , with cos(2Q c ■ Rij)± 1, it is directly 
related with the amplitudes of the exchange modulation. 
Each second-order term that contributes to 7-„-, j.-, also 
adds to 8 Jij, this is explicit in Eq. (|A8|) . As before, 
there are two contributions to jij,]^', the one that sur- 
vives for J^ = J(rij) does not affect the bonds that are 
perpendicular to Q c , but is now symmetric with respect 

toQ c • Q- 

Accounting for the corrections to the bond strength, 
is straightforward and does not require any addi- 
tional consideration. It simply amounts to a change in 
the Fourier-transformed exchange coupling, Jq, which 
determines the ground-state energy and the spin-wave 
spectrum of the Hamiltonian ^,20,21,22,23,24,25,26 The 

significance of this correction, however, is in that a 
change in Jq applies directly to the ground state en- 
ergy, which, therefore, is corrected in the same order, 
0(e 2 ). As we shall see below, the first-order contribution 
to the exchange modulation in Eq. O(e), corrects 

the ground state energy of the spin Hamiltonian only 
in the second order of perturbation. Therefore, except 
for special cases^S these two contributions have to be 
treated equally. 

To summarize, a superlattice distortion @ leads to 
a harmonic modulation of the exchange coupling, with 
either the same wavevector Q c , if it appears as a first- 
order correction to Jy, or with the wavevector 2Q C , if it 
appears in the second order, ~ e 2 . There is also a second- 
order correction to the bond strength. In the most gen- 
eral case these are described by Eq. ®. In what follows, 
we study the effect of the exchange modulation on the 
ground state of the exchange part of the spin Hamilto- 
nian Q (ie the case D = H = 0), which now reads, 



(5) 

Here we have introduced a complex = j^ + ij^, and 

omitted the tildes, keeping implicit that Q c , Jy, j^, j^ 
are all appropriately chosen in accordance with the sit- 
uation, as discussed above. While in the absence of the 
distortion would satisfy all symmetries of the lattice, 
exchange constants in Eq. (JSJ possess only those sym- 
metries of the un-distorted lattice which preserve Q c and 
the polarizations e^ei, 1 (this includes all translations). 

III. CLASSICAL GROUND STATE OF THE 
MODULATED-EXCHANGE HAMILTONIAN 

To proceed with finding the ground state of the 
modulated-exchange Hamiltonian (jSJ), we first rewrite it 
in terms of the lattice Fourier-transforms, taking advan- 
tage of the lattice translational symmetry, 

]y = ^2{ J Q S Q ■ S -q+ 3q+±Q a Sq ' S -q-Q c 

q 

+ jq-±Q c S V S -q+Qj ■ ( 6 ) 

Here j q = J2r tJ Jij exp(-i<7 • r tj ) = j_ q , similar with 
the Fourier-transforms Jq and Sq introduced above. We 
note, however, that unlike Jq which is real, j q is, in 

general, complex, so j„ 7^ j _„■ 

A. General approach 

A general approach to finding the classical ground state 
for a system of N equal spins on a simple Bravais lattice 
that are coupled by isotropic Heisenberg exchange in- 
teraction was developed in Refs. I20l2ll22l and recently 
discussed in Ref. |2(J We need to solve the mathematical 
problem of finding the absolute minimum of a function 
©, which depends on N classical vector- variables, Sq, 
under N constraints, 

s 2 = s 2 yj, (7) 

that are imposed on the length of classical spins. In terms 
of the Fourier-components these constraints become, 

^Sq' • Sq- q > = S 2 5 qfi ,Vq, (8) 

q' 

where 5q,q' is the 3D Kronekker symbol. Upon introduc- 
ing N Lagrange multipliers, Xq, a straightforward varia- 
tion leads to the following equations for the spin structure 
that minimizes the Hamiltonian ©, 

JqSq+ jq-±Q c S q-Q c + 3q+±Q c S q+Q c 

-J2 X 9'Sq-q'=0, (9) 

q' 
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V q, which have to be solved together with Eq. JSJ. 
Substituting these into Eq. ©, we obtain the ground 
state (GS) energy per site, 



which depends only on Xq with q = 0. While this 
suggests searching for the solution with Xq — Ao#g,o, 
it is easily verified by direct substitution that Eq. © 
does not allow any non-trivial solutions of this type for 
jq =/= 0. In fact, using such an "ansatz" for Xq is 
equivalent to replacing (JJJ with a single "weak" condi- 
tion, Ylj &j = NS 2 . This "weak" condition approach, 
also known as the Luttinger-Tisza method, was widely 
employed throughout the early studies of complex spin 
structuresj22i2ii22i£i24 as it does lead to the correct solu- 
tion in several important particular cases, including the 
case of a simple exchange spiral. General reasons for the 
failure of this approach, which, in particular, occurs for 
Hamiltonian (|TJ with D ^ 0, H ^ 0, were discussed in 
Ref. |2y. Current situation, in the form of Eqs. ©, ©, 
presents another example of such failure. 

An alternative approach to finding the ground state 
and the spin-wave spectrum of the Hamiltonian 
for small D and h was used in Refs. It is 

based on a perturbative solution for the real-space spin 
structure that is stationary with respect to small de- 
viations in the form of a slightly distorted flat spi- 
ral. The correction to the spiral winding angle, 9j, 
was obtained in the form of the expansion, S6j = 
J2 n ( a ™ cos ( n< 3 ' r j) + Pn sm ( n Q ' r j)) ! where the coef- 
ficients a n and [3 n are of the order 0((^j-)^ , (7)")- As 
a result, additional harmonics in the Fourier-transform 
of the spin structure, at wavevectors ±nQ, n — 2,3, 
and a magnetization, So, appear. In fact, the same re- 
sult can be obtained from the conditional minimization 
of the Hamiltonian Q outlined above, if a perturbative 
solution is searched in the form of a harmonic expansion, 
x q = E„A„5q in Q, where A„^ ~ 0((^)*,§) ■ \\ n \- U 

and S q = Y,n S nQ S q,nQ> \ S nQ\ ~ °( X \n\-l)- lt is this 
approach, which is both natural, and straightforward to 
apply to the Hamiltonian © in order to obtain the spin- 
wave expansion upon expressing spins through magnon 
creation/annihilation operators^ that we shall employ 
here. 

B. Harmonic expansion for modulated exchange 

Because the modulated-exchange terms in the spin 
Hamiltonian allow the Umklapp processes which couple 
Sq and Sq±Q an d, consequently, couple these Fourier- 
components in the Eq. @, we search for the solution of 
the Eqs. I® in the form of the expansion, 

S Q = H { S Q+nQ c S q,Q+nQ c + S *Q+nQ r 5 q,Q+nQ c } . 

n 

(11) 



where Sq +u q ~ 0(e n ). Substituting this into Eq. © 
it is easy to see that a non-trivial solution requires that 
Xq = unless q ± (Q + nQ c ) — ±(Q + n'Q c ) is satisfied 
for some n and nl . Therefore, the general solution for Xq 
can be written as, 

X Q = J2 XnS q,nQ c > (I 2 ) 

n 

where A n ~ 0(e'™'). Here and below, if the limits are 
not specified explicitly, it is implied that the summation 
extends over all integers (this, in particular, implies tak- 
ing the thermodynamic limit, N — > oo)i22 With (|1 1|1 and 
(|12fl the energy minimum conditions of Eqs. © become, 

J Q+nQ c S Q+nQ c + JQ+(„-1)Q c -Sq + („_i)Q o + 
3*Q+(n+^)Q c S Q+(n+i)Q c ~ Xn ' S Q+(n-n')Q c = °> 

n' 

(13) 

V n. Subsequently, upon substituting (jl 1|) into Eq. JBJj 
the equal-spin constraint is rewritten as, 

XI ( S Q+n'Q c ■ S Q+(n-n')Q c ) = °' (I 4 ) 

n' 

2 Yl ( S Q + (n+n')Q c ■ S h+n'Q c ) = S 2 S n fl, (15) 
n' 

V n. At this point, H13 p -I|15 |l is still a complicated non- 
linear system of equations, and remains such even if we 
retain only the terms ~ 0(e) which determine the lowest- 
order corrections to the simple exchange spiral. 

C. Exchange symmetry 

Further progress into finding the perturbative solu- 
tion to Eqs. (|13fl - i|15|) which would describe a weakly 
distorted exchange spiral, {Sq +71 q ;A n }, is made by 
employing a powerful "exchange symmetry" argument, 
which relates back to the Landau theory of phase tran- 
sitions. It was developed in Ref. |3^ as a basis for the 
unified Lagrangian description of the long- wavelength, 
macroscopic dynamics of the complicated spin systems 
with complex order parameters, including spin glasses. 
Subsequently, this approach was used with great suc- 
cess, in particular, for calculating the low-energy spin dy- 
namics in a variety of situations encountered in the non- 
collinear ground states of CsNiCi3-type triangular-lattice 
antiferromagnetsA2£i2£iSi It is based on a very simple 
observation, that a macroscopic Lagrangian (or a Hamil- 
tonian) of a spin system in a state which is described 
by an order parameter (at T « 0), when expressed in 
terms of the canonical variables that parameterize the 
long- wavelength dynamics of this order parameter, has 
to satisfy all remaining symmetries of the ground state 
(order parameter). Practically, this works as follows. In 
exchange approximation, possible ground states are few, 
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easily classified, and Lagrangians are relatively easy to 
write. Perturbation account for the anisotropy, magnetic 
field, etc, adds terms to the Lagrangian (Hamiltonian) , 
which are expansions in powers of the order parameter, 
and whose general form is essentially determined by the 
above symmetry requirement. 

We extend the exchange symmetry argument to the 
microscopic description of the present paper by noting 
that, so long as the solution of Eqs. l|T5 fl - 1(15 )1 is a weakly 
distorted simple exchange spiral Sj = Sq exp(iQ • rj) + 
Sq exp(— \Q ■ Tj) (ie no independent order parameter 
appears in addition to Sq), and if the perturbation does 
not violate the 0(3) spin symmetry of the initial ex- 
change Hamiltonian, all vectors in spin space, Sq, that 
define the corrections to the initial exchange structure, 
have to be proportional to Sq. In other words, the only 
"selected" directions in spin space which can determine 
direction of spin vectors in the perturbation series are 
those resulting from the spontaneous breaking of the spin 
symmetry that already exists in the non-perturbed sys- 
tem, ie those defined by Sq . This holds for the Hamilto- 
nian J^J), because the modulated-exchange terms preserve 
the 0(3) spin symmetry. Consequently, we write, 

S Q+nQ c = £ nS Q , S Q+nQ = e*_ n S_ Q , (16) 

where e n ~ 0(e n ), and eq = 1. This simplifies equa- 
tions 1|13|) - (|15|) tremendously, as they shall now involve 
only scalar variables e n , A„. In addition, Eq. (|14|l is au- 
tomatically satisfied for (Sq) 2 — 0, which holds in the 
case of a simple exchange spiral. It requires that spin 
vectors of real and imaginary parts of Sq = S' + iS" are 

mutually perpendicular, (S' ■ S") = 0, and have equal 
length. This length is determined from Eq. I|15|). 



\S'\ 



s 2 



along with the following set of conditions on e ? 



(17) 



n'=l 



£n+n'£ n i + £— n /£ — n— n' ) ~ ^t, 



(18) 



which have to be satisfied V n > 0. The energy minimum 
conditions of Eqs. ©, (|13[) become, on account of (|16fl . 



jQ+nQ a e n ^ t A n '£ n _„' 
n'=0 

+ 3Q+(n-i)Q c £ n-l + JQ+(n+i)Q/"+i 

oc 

- ^2 (^n+n>£-n> + A* /£„+„/) = 0, (19) 
n' = l 

V n > 0. Similar equations for n < are solved simulta- 
neously with the above, provided that 



A-|n| — A *n|j 



(20) 



in which case they are simply complex conjugates of 1131) . 

We note, that S-q = S q because Sj are real, so 
Eqs. f2T)| just require that Lagrange multipliers used 
to account for the conditions Q are also real. On the 
other hand, e_ n ^ e*. The solution {e„;A n } of Eqs. 
(JTHJ j US! determines the minimum-energy configuration 
of the equal- length spins through Eqs. (| 1 1 11 - (|lf>|) . (|17|l . 



D. Recursion for the perturbation series and the 
leading-order solution 

Although they superficially look cumbersome, equa- 
tions (fT%|l - (fT§|l are well suited for the perturbation treat- 
ment. Indeed, because 



KS Q+n , Qa ~ O ( e N+l»'l) 



(21) 



it is easy to see that the first line in Eq. (|19|l is ~ e™, 
the second line, except for n = 0, when both terms in it 
are ~ e 2 , contains a term ~ e", and a term ~ e n+2 , and 
the last line, like the last sum in Eq. (|18fl . sums up con- 



tributions 



p n-\-2n' 



n' = 1,2,..., ie overall is 0(e 



n+2\ 



Starting with £o = 1, this defines a set of recursion rela- 
tions which determine {£ n ;A ra } for any given n through 
their values for 1, 2, n — 1. In particular, the leading- 
order, ~ 0(e n ), condition on e n , is, 



-n = - £ (£ n -n>£*- n i) + 0{t 



n+2\ 



(22) 



n' = l 

or, explicitly for the first few orders, 

ei+e*i - 0(e 3 ), 

£2 + £*- 2 = 4 + 0(e 4 ), 

e 3 + £- 3 = £ i( £ 2 - e_ 2 ) + 



(23) 
(24) 

0(e 5 ), (25) 



The leading-order correction to Ao, which determines 
the ground-state energy, appears in the second-order per- 
turbation, ~ e 2 . In the same order appears the intensity 
of the new magnetic Bragg peaks, ~ \Sq±q | 2 , that are 
induced by the exchange modulation. Up to this order, 
we obtain from Eq. I|19|) . 

A o = JQ + [3Q-iQ c - A i) £ -i + (j*Q+iQ c - A i) £ i 
+ 0(e 4 ), (26) 

A i = JQ+iQc + { J Q+Q. - A °) 61 + °( e3 )> ( 2? ) 

which, on account of l|23|) . have the solution, 



£l = 



J Q+Q C + J <3 Q c 2J Q 



+ 0(e d ), (29) 
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with e_i = — e* + 0(e 3 ) (for now, we exclude from con- 
sideration a singular case, Jq+q + '^Q Q = -^q)- The 
ground state energy is then determined by 



the perturbation approach fails extrapolation to Q c = 0, 



where the modulation is absent, and S, 



Q±Q C 



Ao - Jq - 



J Q+Q C 



J, 



Q Q c 



2 J. 



Q 



0(e% (30) 



and, unless Jq_±q — Jq+±Q an d the correction van- 
ishes, it is lower than that of the initial, non-distorted ex- 
change spiral, because Jq is a minimum value of Jq, and, 



therefore, J, 



Q+Q C ^ J QQ C 



2Jq > 0. This is m agree- 



Additional soft regions, such as the lines of soft modes, 
often appear in frustrated spin systems, due to the ac- 
cidental cancellation of the interactions. In such cases 
spin system is extremely sensitive to structural modula- 
tion with the wavevector that is close to the soft region(s) 
of the dispersion. The same is true for distortions that 
propagate along the direction(s) of weak interaction (and 
weak magnon dispersion) in quasi-low-dimensional spin 
systems. 

From Eqs. (|10fl and l|3U|) . the ground state energy is 



ment with a very general argument, that a non- vanishing 
second-order perturbation correction always lowers the 
ground state energy. 



Eqs 
N 



= J Q S 



iQ-hQ, ~3Q+kQ. 



S 4 



-(Qc)VQ 



+ 0(e 4 ). 



E. Some remarks 

It is useful to express the results obtained in the pre- 
vious section in terms of the spin-wave spectrum, 



(35) 



In the general case, it is lowered in response to the ex- 
change modulation. This occurs as a result of the ap- 
propriate adjustment (bunching) of the initial, single-Q 

spiral spin structure through the appearance of the ad- 

u q = SJ2(J q - Jq)(J q+ q + J q q - 2Jq), (31) ditional Fourier-harmonics, Sq +11 q , n = ±1, ±2, . 



and the (/-dependent transverse (perpendicular to the 
spin plane) classical static staggered spin susceptibility, 37 



x±(q) 



1 



2(J q -J Q ) 



(32) 



of the initial, non-distorted, single-Q exchange spiral. 
The leading new Fourier-components of the spin density, 
Sq±q of Eqs. (EH, (EH, are 



= ( 3Q-kQ*~ J QHQA 
' Q+Q < { x±(Qc)" Q S- 2 J 

= ( -^Q-kQc + i*Q+iQc 
' Q ~ Q * \ x±(QcW Qr S- 2 



S Q +0(e 3 ), (33) 



'Q 



0(e 3 ).(34) 



3 


(Q c -V) 


Jq 


2 


(SQ-V) 


(Q c 


V) 2 Jq 


+ 3(Q c -V) 2 Jq 



It is clear now, that the singular case mentioned in the 
previous section corresponds to the exchange modulation 
with the wavevector at which the spin-wave energy van- 
ishes. In this case, unless the numerator in Eqs l|33|) . H34|) 
is also zero, the leading corrections diverge, and the per- 
turbation approach fails. Generally, there are two soft 
spots in the spin-wave spectrum of the exchange spiral. 
They correspond to the Goldstone modes at q = and at 
the magnetic ordering wavevector, q = Q. For Q c = Q, 
the numerator, (J in — j an), is, in general, non-zero 
and the corrections diverge. Unless q = is a special 
(extremum) point of j„, the numerators in Eqs i|33[l and 
(|34fl vanish ~ Q c in the limit Q c — > 0, while uiq in the 

denominator is < Q\. Therefore, for a sufficiently long- 
wavelength distortion the corrections may become arbi- 
trarily large. It is not at all unexpected, though, that 



In 

addition, the pitch of the primary spiral component, Sq, 

may also change, Q — > Q, because the spiral propagation 
vector, Q, is now defined by the minimum of the cor- 
rected energy, Eq. (fSTTf) or ij^ST) . In the case where the 
modulation of the crystal structure is long-periodic, ie 
for small Q c <C 1, and assuming that correction to the 
propagation vector, SQ = Q — Q, is small, \5Q\ <C Q, 
we can expand Jq ± q an( l J q±iq m Taylor series and 
obtain, 

Ao = Jq + \ (SQ ■ V) 2 J Q 



+ 0(e 3 ), 
(36) 



where V = -§g. Note, that unlike Eqs {3"3"J). this, gen- 
erally, does not diverge for Q c — * 0. If, in the correction 
term, we cancel Q\ and then expand the denominator 
keeping only the leading term in SQ, Eq. takes the 
form A = Jq + c + cx(SQ) + c 2 (SQ) 2 , with c x = 0(e 2 ) 
and C2 = O(e). Clearly, minimum of this expression oc- 
curs, in the general case, for non-zero SQ = 0(e). How- 
ever, if the linear term ~ SQ in the denominator of Eq. 
1361 vanishes, then SQ — 0. Also, in many important cases 
such as, for example, the nearest-neighbor non-frustrated 
antifcrromagnet, Q is a special symmetry point of j q 
(Jq ~ cJq), and the correction term identically van- 
ishes by symmetry. Therefore, simple exchange struc- 
tures, such as antiferromagnets, are usually insensitive 
to small long-periodic modulations of the exchange cou- 
pling in the spin Hamiltonian (J5J. 
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IV. SOME EXAMPLES 

Now we shall apply the formalism developed above to 
several representative one- and two-dimensional (ID and 
2D) systems. The fact that the ordered mean-field (MF) 
ground state, such as analyzed in this paper, in low di- 
mension is unstable against the fluctuations^ will not 
be a concern here. First, a single-Q c structural distor- 
tion considered in this paper is homogeneous within the 
crystallographic planes that are perpendicular to Q c . In 
many cases these planes contain one or two unit cell direc- 
tions, so the modulation only exists along the remaining 
direction(s), and the distortion is explicitly 2D, or ID, re- 
spectively (cf, Fig. 1). The bonds that are not changed 
by the distortion cancel out in the resulting expressions 
for the corrections to the spin structure and to the ground 
state energy, Eqs. if^l) . i!5U|l . Therefore, even though the 
Hamiltonian Q may be on a 3D lattice, the distortion 
corrections will be the same as for the lower-dimensional 
system. Also, in many quasi-low-D materials the essen- 
tial physics is ID or 2D, while the MF order is stabilized 
by weak interaction in the remaining direction(s). Such 
is the case of high-T c cuprates which are made of two- 
dimensional square-lattice layers. In the presence of a 
charge order the layers are modulated, as illustrated in 
Fig. 1 (b). In fact, the MF analysis has certain value 
even for purely low-D systems, because it highlights pos- 
sible phases and instabilities, and guides the behavior of 
the critical points, at least asymptotic (S — > oo). 



A. n-merized ID antiferromagnetic chain 

As a simple example, consider the antiferromagnetic 
spin chain with ri-periodic nearest-neighbor exchange 
coupling (ferromagnetic case, J < 0, is trivial, because 
small bond modulation has no effect on the ground state). 
It corresponds to a modulation with Q c = — , Fig. l(al). 
In this case, Jq = 2 J cos q, J„ = 2jcosq, the modified 
GS energy is determined by 

A = 2 J cos Q (l + [7/jf tan 2 q) , (37) 

and does not depend on n. Clearly, there are new local 
extrema at cos 2 Q ~ \j / J\ 2 , corresponding to a spiral 
with pitch of 7r/2, which appear on account of the ex- 
change modulation. However, the new minimum is very 
shallow, Ao ~ — 4|j'|, much higher in energy than the 
global minimum, Ao = — 2J, which remains at Q = it. 

A more interesting ID example is illustrated in Fig. 
I(a2), where in addition to the nearest-neighbor coupling, 
Ji, there is a frustrating next-nearest-neighbor exchange, 
,/ 2 • In this case Jq = 2 Ji cos q + 2 J 2 cos 2q and, for 4J 2 > 
Ji, frustration (for classical spins) results in a spiral MF 
ground state, with a pitch defined from the condition 
cos Q = — jj^ . For the sake of simplicity we consider only 
modulation of the nearest-neighbor coupling, ie j 2 = 0. 



In this case, as before, j„ = 2j 1 coaq, and, 

2 Ij'J 2 sin 2 Q 
~ Jq + J lC osQ + 4J 2 cos 2 ^cos2Q' (38) 

The minimum of this expression is achieved for Q = Q + 
SQ which differs from Q = cos -1 I — ) . For small 8Q 




FIG. 1: Examples of the modulated exchange patterns dis- 
cussed in the text. Different bond thicknesses illustrate dif- 
ferent coupling strengths, (al) n-merized ID chain with mod- 
ulated nearest-neighbor coupling, Q c = — ; (a2) dimerized 
(n=2) antiferromagnetic chain with frustrating next-nearest 
neighbor coupling, 4J2 > Ji. (b) "Stripes" on a square lat- 
tice with diagonal modulation, Q c — (2£ , — ); n = 4 case is 
shown, (c) Generalized "staggered row model" on a triangu- 
lar lattice, obtained with Q c = ( — ,—). Shown is the case of 
n — 4; n equals 8 in the practical case of RbMnBr3. 
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the leading correction is, 



SQ = 



h_ 

2J 2 



Q 



tan Q(l- cos^cos2Q) 



(39) 



and diverges for Q — > 7r, ze for frustrated ID antifcr- 
romagnet with J 2 = 0.25Ji. Perhaps, this interesting 
finding is an indication of an instability towards a spon- 
taneous n-merization and a spin gap formation in the 
vicinity of this point. In fact, for S=l/2 quantum spin 
chain the whole region J 2 > J2 C w 0.25Ji is believed to 
belong to a spin-gap phasef22iiMi including the special 
point, J 2 = 0.5 Ji, where the spontaneously dimerized 
ground state is known exactly^ While the most recent 
numerical estimate for S=l/2 chain is J 2c ~ 0.24Jipi a 
J2c = 0.25 Ji for S — > oo coincides with the earlier result 
obtained from the semi-classical mapping on the nonlin- 
ear sigma-modeliSi Clearly, studying the susceptibility 
of spin system towards bond modulation is a proper way 
to investigate its instability towards n-merization (pla- 
quette formation) and to characterize the corresponding 
phase diagram. To this end it seems possible to develop a 
perturbation approach similar to the one presented here 
and starting with the Hamiltonian (jSJ) also for quantum 
spins, but this goes far beyond the scope of this paper. 



B. Square lattice antiferromagnet with diagonal 
modulation 



A square-lattice, nearest-neighbor antiferromagnet 
with diagonal modulation corresponds to Q c = (—,—), 
Fig. 1(b), and may be of direct relevance for the 
charge-ordered stripe phases in doped LSCO cuprates 
and related perovskites. In the isotropic case the bond 
strengths and the modulation amplitudes are equal in 
two directions, Jq = 2 J (cos (q ■ a{) + cos (q ■ a 2 )), j q — 

eJq, and, upon switching to Q = q ■ ( ai + a2 ) anc j Q> = 
q ■ ( ai ~ Q2 ) j the problem is factorized and reduced to the 
one-dimensional one considered above. The distortion- 
corrected spin GS energy is determined by Ao of Eq. 
(|37|l . which is simply multiplied by cosQ'. As before, 
we find that nearest-neighbor antiferromagnetism is sta- 
ble with respect to the bond modulation. In fact, the 
same conclusion is reached even if the amplitude of the 
bond modulation is different in two directions, so that 

jq = 2 Ji cos (<? ' a i) + 2 i 2 cos (l ■ fl 2)- 

Again, an interesting situation occurs if there is frus- 
tration. In the case of the square-lattice antiferromagnet 
it is introduced by the diagonal, next-nearest-neighbor 
coupling, J' > 0. This model has been extensively stud- 



ied in recent years t 



43.44.45.46.47.48 



since it was predicted 



that in a region of parameters in the vicinity of J/ (2 J') 
I it has a disordered, spin-liquid ground statei 43 ' 44 For 
S=l/2 quantum spins this was proposed as a possible 
candidate for the resonating-valence-bond statei 4 ^* 5 ^ In 
this case, Jq = AJcosQcosQ' + 2J'(cos2Q + cos 20/), 
and, if both side and diagonal bonds are modulated, 



jq = 4j cos Q cos Q' + 2j (cos2Q + cos2Q'). Conse- 
quently, we obtain from Eq. (|3(J|) . 



4sin 2 Q 



j cos Q' + 2 j cos £ cos Q 



A ° Jq+ JcosQcosQ' + 2J'cos 2 ^cos2Q ' ^ 

In the absence of bond modulation the ground state 
is determined by the hierarchy of the local minima of 
Jq. It depends on the relative strength of the nearest- 
neighbor coupling J, and the next-nearest-neighbor, di- 
agonal coupling J', which is parameterized by a = -^jj. 
For weak frustration, a > 1, the global minimum is that 
with sinQ = sinQ' = 0. It corresponds to the conven- 
tional, collinear Neel antiferromagnetic order with a sin- 
gle propagation vector Q = (tt, tt), and the ground state 
energy is j^E^^) = — 4JS' 2 (1 — ^). Although there are 
four equivalent Q-points in the BZ, (±7r, ±7r), (^-^,±71"), 
which restore the lattice C4 rotational symmetry, they 
are related through addition of the appropriate recipro- 
cal lattice vectors r, so there is no true GS degeneracy in 
the Q-space. The only GS degeneracy is the rotational 
symmetry in spin space which corresponds to the 0(3) 
symmetry of the Heisenberg spin Hamiltonian. 

For strong frustration, a < 1, additional ground state 
degeneracy occurs on the MF level. For a < 1 there are 
two non-equivalent lowest-energy minima of Jq , satisfy- 
ing cosQ = cosQ' = 0. They correspond to two pairs 
of equivalent Q-points in the BZ, (±7r, 0) and (0,±tt), 
which represent the antiferromagnetic order propagat- 
ing along x and y axis, respectively. The GS energy is 



-AJS 2 ^-. This double degeneracy 



in Q-space can be used to construct a continuum of states 
which are the linear combinations of the above two. This 
continuous GS degeneracy is usually described in terms 
of two decoupled y2 x \[2 sublattices based on the diago- 
nals of the original square lattice, which is transparent in 
the limit J' ^S> J. Each sublattice has an antiferromag- 
netic order, but there may be an arbitrary angle between 
the two, because the mean field from one sublattice can- 
cels on the sites of the other. This continuous degeneracy 
is lifted by zero-point or thermal fluctuations which pre- 
fer collinear arrangements of the two sublattices in the 
GS. This is a famous example of the order by disorder 
phenomenon in a frustrated magnet. 51,52 

The most interesting situation occurs for a = 1 , when, 
on the MF level, there is a continuous GS degeneracy 
even in the Q-space. The minimum condition for Jq 
becomes cosQ = cosQ', and is satisfied for any spiral 
with the propagation vector Q that belongs to the square 
with the vertices at (±7r, ±7t), (T 71 ": ±7r). All these states 
belong to the global minimum, and have the same energy, 
±E a=l = -2JS 2 = -4J'S 2 . It is this continuous Q- 
space degeneracy which is at the origin of the spin-liquid 
phase conjectured for the values of a in a finite region 
around the special point a = 1 

Importantly, for a < 1 the spiral states with Q w (jr, w) 
are in close competition with the collinear states. In par- 
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ticular, consider an extremum of Jq which is a local min- 
imum along the diagonal direction, parallel to the lattice 
modulation wavevector Q c , whose energy in the absence 
of the modulation is jjEq = —2aJS 2 . It corresponds to 
a spiral with the propagation vector defined by Q' = 0, 
cosQ = —277, ie Q = (cos~ 1 (— ■%jr), cos _1 (— ijr)) (there 
is also a degenerate state with Q at 90°, respecting the 
Q <-» Q' symmetry of the square lattice). Except for 
a = 1 the energy of this extremum is higher than that 
for the decoupled antiferromagnetic sublattices, Er^m. 
It is clear, however, from Eq. (|4U|I that, while the energy 
of the collinear antiferromagnetic states is insensitive to 
bond modulation, the energy of the spiral GS can be low- 
ered as it adapts to distortion! Therefore, at least on the 
MF level, spiral may become the lowest energy state (ie 
the ground state) for some range of the parameter a in 
the vicinity of 1 (whose width is ~ 0(e 2 )). For the long- 
periodic modulations, Q c -C 1, and for j = 0, it is easy 
to find that spiral phase is stable for 1 — \j /J\ 2 < a < 1. 
The "nominal" spiral propagation vector Q is obtained 
by minimizing Eq. H4U|) . similarly with the case of the 
frustrated ID chain. Again, our finding clearly indicates 
the instability of the frustrated square-lattice antiferro- 
magnet with J j (2 J') close to 1 with respect to the bond- 
modulated states. By selecting spiral spin GS the "or- 
der by distortion" mechanism proposed here competes 
with the "order by disorder" phenomenon, which prefers 
collinear states. At least for large enough S spiral always 
wins in some vicinity of a = 1 . 



shown in Fig. 1(c), and the one with n = 8 is relevant 
for the phases realized in RbMnBr 3 and KNiCl 3 i 12 i 13 i 14 
In this case j„ = Aj cos Q cos Q' + 2j cos 2Q and it is easy 
to see that the energy of the modulated state is given by 
the same expression as for the frustrated square lattice, 
Eq. H40fl . but with j — j and J' — J (as in the previous 
example, we use Q — q ■ ( Ql + a2 ) ail d Q' = q ■ ( Ql ~ a2 ) ), 
Exchange modulation leads to a deviation from the com- 
mensurate 120° triangular spin structure. The leading 
correction to the propagation vector is determined from 



cos Q 



+ S, where 



cos — 



cos^) 



In the two limiting cases, 
2 



■ cos - 



§ cos 2 



-) 



6 = 



and 5 « (£) 



(41) 

2 and n 3> 1, we find 
respectively. It is clear 



that S is small for all n, so there is no evidence for an 
instability towards n-merization in the case of ideal tri- 
angular lattice. Perhaps, such evidence can be found in 
the anisotropic, quasi-lD case, where exchange in zigzag 
rows is much smaller than that in straight rows, or vice 
versa. 



SUMMARY AND CONCLUSIONS 



C. Generalized row models on triangular lattice 

Because of inherent frustration, the triangular lat- 
tice antiferromagnet, Fig. 1 (c), is a very interesting 
case to consider. Without modulation, a minimum of 
Jq = 2J(cos (q ■ a±) + COS (q ■ a 2 ) + COS (q ■ (di + a 2 ))) 
is achieved for a commensurate spiral GS with propa- 
gation vector Q = (^f,^f)- It corresponds to a non- 
collinear GS structure where spins are aligned along one 
of the 3 directions, at 120° degrees with each other. 
Structural distortion may result in a variety of coupling 
patterns where equivalent bonds are related by transla- 
tions perpendicular to Q c . These can be classified as 
"generalized row models", the simplest of which is the 
original row model of Ref. @. It can arise, for exam- 
ple, as a second-order effect from the distortion with 
Q c = (±7r,7r). Because 2Q C = r (t is a reciprocal lat- 
tice vector), it is equivalent to a homogeneous uniaxial 
compression, ie it simply changes Jq and does not result 
in a bond modulation. Nevertheless, the ground state 
becomes incommensurate, with pitch of the spiral deter- 
mined by the coupling anisotropyi&£ 

A general modulation of the triangular lattice with 
Q c = (^jS - ^f) results in a "staggered row model", 
where the horizontal rows have equal coupling. On the 
other hand, modulation with Q c = ( — , — ) leads to a 
"zigzag row model" , an example of which with n = 4 is 



In a great variety of important practical cases the com- 
plex crystal structure which is at the origin of the in- 
tricate magnetic behavior in magnetic material results 
from a small superstructural distortion of a much sim- 
pler structure, in which the magnetic ions form a prim- 
itive Bravais lattice. Reduction of the crystal symmetry 
related to the appearance of even a single, commensurate 
with the original lattice, superstructural Bragg reflection 
at a wavevector Q cl formally requires folding the original 
Brillouin zone to a much smaller one, many times re- 
ducing its volume. Consequently, the magnetic system is 
usually described in terms of multiple spin sites, and mul- 
tiple sheets of spin excitations. Not only does this greatly 
complicate understanding and predicting magnetic prop- 
erties, such violent modification of the q-space clearly 
seems an unsatisfactory way to account for a small dis- 
tortion of the crystal structure. Moreover, BZ folding is 
not an option for the incommensurate structural modu- 
lations, such as arise in various charge-density-wave or- 
dered states. 

In this paper the effect of a small lattice modulation 
with single propagation vector Q c on the system of local- 
ized spins, coupled by Heisenberg exchange interaction 
||TJ, was considered. It was found that lattice distortion 
results in a modulation of the exchange coupling which, 
to the leading order, is parameterized by no more than 
two constants per bond, Eq. JSJ (this is valid irrespective 
of whether the spins are quantum or classical) . There are 
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also corrections of the order 0(e 2 ) (e is a small parameter 
that parameterizes the lattice distortion) to the exchange 
constants in the covariant part of the Hamiltonian. It 
should be noted here that, although the distortion con- 
sidered is small, the resulting corrections to the coupling 
constants need not be small compared to the initial val- 
ues of the couplings in the spin Hamiltonian, which may 
be small themselves. Therefore, the result expressed by 
Eq. (J3J) is quite general and does not automatically imply 
the condition <C |Jy|. In other words, the "Umk- 
lapp" terms ~ SqSq±Q added to the spin Hamiltonian 
by small lattice distortion may be relatively large, even 
larger than the original exchange interactions. 

While it would be interesting to study the modulated- 
exchange Hamiltonian for quantum spins and for the ar- 
bitrary values of ■ • / Jy | , it is a formidable task which is 
beyond the scope of this paper. Here we developed a per- 
turbative scheme for finding the mean field ground state 
of the Hamiltonian (JSJ) which is valid for classical spins, 
S 3> 1, and in the case of small exchange modulation, 
\jij/Jij\ ~ e < 1. One of the initial motivations here 
was to develop a basis for the spin- wave theory in weakly 
distorted crystal structures in terms of the modification 
of the magnon spectrum in the original, large Brillouin 
zone of the non-distorted Bravais lattice. The other, no 
less important motivation, was to see whether it would 
be possible to understand, already on the mean field 
level, the incommensurate phases observed in the dis- 
torted triangular-lattice antiferromagnets RbMnBr3 and 
KNiCi3, and in the doped, distorted square-lattice anti- 
ferromagnets, such as LSCO or related Ni, Mn and Co 
materials, such as Lai. 5 Sr .5CoO4iS 

The essential results of this paper are expressed by 
Eqs. HI}, (©, or, equivalently, by Eqs. - (|55)> . 

They show, that a transverse, equal-spin spiral structure, 
which is the ground state of the initial Heisenberg Hamil- 
tonian, adapts to the exchange modulation through ap- 
pearance of the additional Fourier-harmonics, Sq +u q , 
n = ±1, ±2, ... (bunching). As a result, in a general case 
the GS energy is lowered by exchange modulation. In ad- 
dition, the pitch of the primary spiral component, Sq, 

may also change, Q — > Q, because it is now defined by 
the minimum of modulation-corrected energy, Eq. 1|35[) . 

Applying these results to several particular examples 
of the topical frustrated spin systems appears quite re- 
vealing. We find that in the case of the frustrated 
square-lattice antiferromagnet with diagonal coupling J', 
such that a = J/(2J') is close to 1, lattice modulation 
opens a region of stability of the incommensurate spiral 
phase. This "order by distortion" phenomenon^ com- 
petes with "order by disorder" , which prefers collinear 
arrangements of two antiferromagnetic sublattices. The 
incommensurate spiral phase with the propagation vector 
Q = (n±5, it ±5) close to (it, n) wins for the range 0(e 2 ) 
of the parameter a in the vicinity of a = 1. This provides 
a plausible explanation for the incommensurate spin or- 
der observed in Lai.sSro.sCoO^^ and in a number of 



other doped perovskites, and may also be of a direct rel- 
evance for the doped LSCO materials. It is important to 
mention here that incommensurate spin-ordered phases 
are among the most interesting and puzzling features of 
doped layered perovskites. In the absence of distortion 
one needs at least a third-neighbor coupling in order to 
stabilize spiral MF ground state for the Heisenberg spin 
Hamiltonian on square lattice. 

Exchange modulation in Heisenberg antiferromagnet 
on a distorted triangular lattice leads to an incommensu- 
rate shift of the spiral propagation vector, in qualitative 
agreement with what is observed in RbMnBr3. How- 
ever, Eq. (|41|) implies that Q is decreased compared 
to Q = (^ , 2 3 E ) for the ideal triangular lattice, while 
Q sa 2ir ■ (0.357, 0.357) is observed in experimentsi l: ?i 14 
Therefore, it is likely that the shift in RbMnBr3 is mainly 
due to the anisotropic corrections to the nearest-neighbor 
coupling, which are captured already in the simplest row 
model. 8 ^ Nevertheless, correction of Eq. I|41[l is not 
unimportant. In the case of n — 4, which may be rel- 
evant for RbMnBr3, it gives the same magnitude shift of 
the ordering wavevector, SQ, as measured in experiment, 
for \j/J\ 2 « 0.2 (for n = 8, | j/J\ 2 « 1 is needed). 

More importantly, bunching of the spin spiral as a 
result of the lattice distortion provides, already on the 
mean field level, a possible explanation for the com- 
mensuration transition in RbMnBr3 and for the long- 
periodic lattice-commensurate structure in the related 
phase of KN1CI3. Indeed, an easy-axis anisotropy, such 
as in KNiCl3 , or a magnetic field applied within the easy 
plane, as in RbMnBr3, also lead to bunching of the ex- 
change spiral, generating additional Fourier- harmonics, 
Sq, at q = nQ, n — ±2,±3,.... 26 Appearance of these 
Fourier-components in the spin distribution lowers the 
spin-anisotropy and the Zeeman energy, but competes 
with the modulated exchange, which requires additional 
Fourier-components at q = Q + nQ c . Therefore, for some 
finite value of the easy-axis anisotropy, or the in-plane 
magnetic field, a commensuration transition may be ex- 
pected, where Q becomes equal to mQ c with some inte- 
ger m. In the lattice-commensurate phase both sets of 
additional harmonics coincide, and both the modulated 
exchange energy, and the spin-anisotropy and Zeeman 
energy, can be lowered simultaneously (this, of course, 
should off-set the increase in the un-modulated exchange 
energy caused by the shift in Q). Extending the results 
of this paper and those of Ref. |26| to the Hamiltonian 
with modulated exchange of the Eq. (0 and D, H ^ in 
order to map out such phase diagram is one of the most 
obvious directions for further studies. 
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APPENDIX A: EFFECT OF LATTICE 
MODULATION ON EXCHANGE COUPLING 

The simplest assumption, and the one which is most 
often employed in literature^ is that the exchange cou- 
pling Jij between the magnetic ions (spins) at positions 
r, and rj only depends on the distance between the sites, 
Jij = J{\ r ij\)- However, in magnetic dielectrics this cou- 
pling most often results from the superexchange and, 
therefore, generally also depends on the position(s) of 
ligand ions which bridge the superexchange path, 

Jij = J(\r^l{r^}) = J{n 3l {u^}). (Al) 

Here fXij numbers the ligands that participate in the 
ij bond and r MiJ are their positions which are most 
naturally parameterized in terms of the offsets, u^ iJ = 



Rij, from the bond center, R^ — ^(Ri 



R-j 



In addition, the superexchange coupling may also de- 
pend on the angles of orbital overlaps. This dependence 
can be parameterized in terms of the positions of some 
particular symmetry points in the local electron density 
distribution, and accounted for in IjAljl by including these 
points among {r^ ij }. While these additional degrees of 
freedom do lift some non-essential symmetries which are 
present in the particular case when the exchange cou- 
pling only depends on the bond length, = J(\rij\), 
they do not change the general structure of the correc- 
tions to the exchange coupling resulting from the lattice 
modulation which are summarized by Eq. @ in the main 
text. In what follows, we first discuss the particular case 
of Jij — J(\rij\), and then the general case of Eq. IjAljl . 



drij 

is given by the first-order term, 



For (eij-rij ) ^ and g ^ r ' ] ' ^ the leading contribution 

e, and we obtain, 



(A2) 

where, as usual, = Tj — n and Rij = \{Ri + Rj)- 

Expanding the exchange coupling, = J(r^), mod- 
ified by the distortion (jA2j) in a Taylor series in small 

Ti, 



displacement, e,j <C iy, 



^■^—) Jinj), (A3) 



dr 



ij 



we find, up to a second order in (^^j ~ K< 1, 



J' -J 



><.,, ■ r,,i <)Jir :J ) 1 (eij ■ r,j) 2 d 2 J{r i3 ) 



dry 2 r 2 - Drfj 



2 I ra rf- j <9r„ 



+ 0{e 3 ). (A4) 



J'in = J *3 + hi C0S (Qc ■ Rij) - j'L sin(Q c • R t 



0(e 2 



(A5) 



where 



f ^(Ss^^M^.r.j). ,A(i. 



J: 



-2 sin 



2 / ^ij ^^ij 

^lM.\^pl {ei . rij) . (A7) 



Clearly, to the first order, there is no Q c -independent 
correction which would change the bond strength J^. 
Such a correction does appear in the second order in e. 

If (e.y • rij) — 0, ie the displacements are perpendic- 
ular to the bonds, the leading correction is given by the 



second-order term 



In fact, this important situa- 



tion is often encountered in practice, in particular, it is 
the basis for the so-called "row models" for triangular 
lattice antiferromagnets. These are typically thought to 
be realized through a modulation of the hexagonal lattice 
where the displacements of the magnetic ions are parallel 
to the Cq symmetry axis, and are perpendicular to the 
bonds in the hexagonal plane in which the modulation 
propagates. RbMnBr3 is believed to present an example 
of such situation. For such case of a transverse structural 
modulation we obtain, 

J'ij =Jij + Jij cos(2Q c • I^j) - j", sin(2Q c • Ry) 



+Sij [ {hi? + ) 2 ) 2 + • rv,) 2 , e 3 ), (A8) 



1. Modulation of the bond length only 

First, consider the effect of displacement of the mag- 
netic ions alone. In presence of the superstructural mod- 
ulation J21, the bond lengths become = \nj + £y|, 
with eij — Ej — 6j given by 

€;■, = -'sin [ c J v ) (-ei sin(Q c • Rij) + e 2 cos(Q c • R {j )) , 



'ij — U V ^ Jij ,MO \*'*c *HjJ Jij °"H**tfc 
{'Jij 

whore Sij = sign ( fg^f p and 



J 



= sm 



2{Qc-r ij \dJ(r ij ) ( £ 2_ e 2^ (A9) 



1 — sm 



Tij dfij 

•„2 ( Qc ■ rij_\ dJjjTjj) 



2(ei-e 2 ). 



(A10) 



The first-order correction to J[j also vanish if aJ J^ 13 - = 0. 
In this case the leading contribution comes from the last 
second-order term in Eq. ljA4jl . and is also expressed by 

Eq. (|A8j) . but with = sign ( ^P^ j, and 



'" " - 2 l ' Qc ' rij ^ d2Jinj) ((e 2 -rij) 2 - (et-njf) , 



3 ti =sm 



2 J ^drl 



(All) 



f!j^(9^M.\ £^ 2 , €| . r .. )I63 . r; ,, 



(A12) 



in place of ljA9H . (jA10jl . 

A Q c -independent correction which changes the bond 
strength J^ first appears in the second order, Eq. 
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(|A8|) . In general, it is obtained by summing up all Q c - 
independent contributions from all second-order terms 

in l|A4|) . Their common multiplier, sin 2 ^ , makes 

the structure of this correction rather simple. It does 
not affect the bonds that are perpendicular to the direc- 
tion of propagation of the lattice distortion, while those 
bonds that are symmetric with respect to this direction 
are modified equally. On the other hand, j jt and jU of 
Eqs. IjAfijl . (|A7|1 which describe the first-order correc- 
tion to the exchange coupling resulting from the lattice 

modulation on are ~ sin ( -^if 12 ^ , and, therefore, are 
antisymmetric with respect to Q — > — Q. 

2. Modulation of the ligand positions, etc 

It is also straightforward to account for the dependence 
of the superexchange coupling, Eq. IjAlfl . on the posi- 
tions, , of the ligand ions and the symmetry points 
of the local electron density distribution which define the 
orbital overlaps. In the presence of the superstructural 
distortion JSJl, 

(r^i)' = r n*i + e £cos(Q c • r***) + sin(Q c • r"«), 

(A13) 

where fi indexes different types of r Mij positions within 
the unit cell, and the polarization vectors e^ 2 which pa- 
rameterize the displacement for point of type \x are deter- 
mined by the superlattice Bragg intensities that appear 
with distortion. This can be rewritten as, 

(u^)' = + ef 3 cos(Q c • Rij) + e% ij sm{Q c ■ Rij), 

(A14) 

where the new polarization vectors, *| now depend on 
Q c . They are obtained by rotating e^ 2 through an angle 
4>fj,ij — (Q c ' u ^ ij )i and subtracting e lj2 cos(iQ c -ry) (this 
accounts for change in the bond center position Rij), 
correspondingly. 

Consequently, in the general case of Eq. l|A"Tjl . the 
Taylor series (IA3(I for JL has to be amended, by adding, 



n fj.ij * ' 

. (A15) 

whose first and second-order terms are easily rewritten in 
the form of Eq. i|A5|l and Eq. (|A8|) , respectively. There- 
fore, the account for the modulation of the positions u^ i] 
in the general expression for the superexchange, 1A1|I . 
simply amounts to amending the coefficients and jU 
in equations (|A5|) and l)A8|l . in accordance with Eqs. 
I|A14|I . I)A15(I . For example, additional first-order terms 
which appear in (|A15|) change the expressions of Eqs. 
(|A"6)l and l|A7j) as follows, 

f, • jS~) j ir ^ u " l; ' (A16) 

fj-ij 

hi >f V E • 9^~) J(r«,{«^}) • (A17) 

Using Eqs. (|A14|) . l)A15jl it is easy to write out similar 
expressions for the coefficients of the second-order con- 
tribution of Eq. HA8fl . 

Clearly, a number of symmetry properties of the co- 
efficients jU and jU given by Eqs. (|A6f> - (|A12fl which 
are present for Jjj = J {rij), disappear upon account for 
the modulation of the positions u^'i . In particular, for 
the first-order corrections to vanish, not only should the 
displacements of the lattice sites (magnetic ions) be per- 
pendicular to the bonds, but all of the displacements e^'l 
should be perpendicular to the corresponding gradients 
of Jij — J(r.y , }) with respect to However, 
it is clear that, because jU, jU and, in general, Jy, are 
functions on the non-distorted lattice which also depend 
on the modulation wavevector Q c and the polarizations 

2 , they are invariant with respect to all symmetry op- 
erations of that initial lattice which do not change Q c and 

2 ■ Importantly, this includes the translation group of 
the non-distorted lattice, which means that the new cou- 
plings and the exchange modulation amplitudes which 
are related by that lattice translations are equal. 
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